Theory of spin-conserving excitation of the N — V center in diamond 
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The negatively charged nitrogen-vacancy defect (N — V~ center) in diamond is an important 
atomic-scale structure that can be used as a qubit in quantum computing and as a marker in 
biomedical applications. Its usefulness relies on the ability to optically excite electrons between 
well-defined gap states, which requires clear and detailed understanding of the relevant states and 
excitation processes. Here we show that by using hybrid density-functional-theory calculations in a 
large supercell we can reproduce the zero-phonon line and the Stokes and anti-Stokes shifts, yielding 
a complete picture of the spin-conserving excitation of this defect. 
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Quantum computing and its many exciting applica- 
tions relies on the successful realization of quantum logic 
bits (qubits) that can operate under practically feasible 
conditions. Few physical systems can meet the require- 
ments of controlled quantum coherence and robust op- 
erational conditions. One of the most promising candi- 
dates is the negatively charged nitrogen-vacancy defect 
(N - V~ center) in bulk diamond [J, 0]: the spin state 
of this defect can be manipulated using excitation from 
the 3 A 2 ground state to the 3 E excited state by optical 
absorption (Fig. |T|) . The main advantage of the N — V~ 
center is that it can operate at room temperature as a 
solid state qubit 0, 0, H, @, 0, 0, @|. Detailed under- 
standing of this excitation process is crucial in the re- 
alization of qubits based on diamond. However, achiev- 
ing this level of understanding stretches the capabilities 
of theoretical methods that are usually applied to the 
study of defects in solids because of the special nature 
of the N — V~ defect. Specifically, this defect combines 
strong coupling of ionic and electronic degrees of freedom 
with a many-body character of the electronic states, and 
its interpretation is further complicated by contradictory 
experimental measurements of the emission spectrum at 
different temperatures filficij]. 

In this Letter, we report a theoretical investigation 
of the radiative transitions of the N — V~ center us- 
ing electronic structure calculations, which give a con- 
sistent and accurate account of the excitations observed 
and provide a plausible resolution of the experimental 
situation. We utilized the HSE06 screened Hartrce Fock 



hybrid exchange-correlation density functional [111, [12 1 
to determine the geometry and excitation energies of the 
N — V~ center, and we compare the results to the tradi- 
tional PBE 13| exchange-correlation density functional 
and the experimental data. We find that - in contrast 
to PBE - the HSE06 functional (which reproduces the 
band gap of diamond within 0.5%), can also reproduce 



both the zero-phonon line and the Stokes-shift quanti- 
tatively (within 1.5% in this case). This result demon- 
strates that hybrid functionals improve not only the ex- 
citation of the extended system (gap) but also localized 
ones. This promises a very significant advantage in de- 
fect calculations. Motivated by the success of the hybrid 
functional to reproduce key experimental measurements 
of the excitation process, we calculate an anti-Stokes shift 
of 0.217 eV, measured to be 0.185 eV at usual experimen- 
tal conditions at low temperature From this result, 
we argue that the anti-Stokes shift of 0.065 eV measured 
at low temperature under high-energy-density laser illu- 
mination [lfjj, is most likely due to the local heating of 
the sample caused by the focused laser beam. 

Previous density functional theory (DFT) calculations 
have shown 14, 15, 16|, 17 1 that well defined defect lev- 
els appear in the band gap due to this defect: a fully 
occupied a± level and a doubly occupied two-fold degen- 
erate e-level at a higher energy. The electrons have par- 
allel spins on the e-level preserving the C3i, symmetry of 
the defect; a many-body wavefunction built from single- 
particle states in the gap represents the 3 ' A 2 ground state. 
The excitation can be understood as promoting one elec- 
tron from the a\ level to the e level resulting in a new 
many-body excited state, 3 E, as shown schematically in 
Fig.ffl 

Excitation changes not only the electron wavefunction 
but the atomic structure of the defect as well. Hence, 
the ground state and the excited state will possess dif- 
ferent potential energy surfaces (PES) and different vi- 
brational states, as shown schematically in Fig. [2l The 
transition between the lowest PES will result in the zero- 
phonon line (ZPL) both in absorption and emission, a 
process in which no real phonons are involved in the ex- 
citation or de-excitation process. The ZPL was measured 
at 1.945 eV (yellow light) both in low temperature ab- 
sorption and emission At liquid-nitrogen tempera- 
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FIG. 1: (Color online) (a) The structure of N — V~ center 
in diamond; only first- and second-neighbor C (cyan spheres) 
and N (blue sphere) atoms to the vacant site are shown. The 
yellow and red lobes are contours of the calculated difference 
in spin density for the 3 E state as obtained by the HSE06 and 
PBE functionals. (b) Schematic diagram of the defect states 
in the gap and their occupation in the 3 A2 (ground) and 3 E 
(excited) states. 



ture a broad phonon side band was measured in absorp- 
tion with phonon-related peaks at approximately 2.020, 
2.110, 2.180 and 2.250 eV, with the highest intensity at 
2.180 eV (green light), while in the emission band the 
first phonon sidebands, better resolved than in absorp- 
tion, are found at 1.880, 1.820, and 1.760 eV, with the 
1.760 eV peak (red light) having the highest intensity Q. 

The Franck-Condon approximation is commonly used 
to interpret the excitation spectrum, that is, assuming 
that the electronic transition is very fast compared with 
the motion of nuclei in the lattice. In addition to the 
Franck-Condon assumption, three other approximations 
are commonly assumed (see Fig. [2]). The first is that each 
lattice vibrational mode is well described by a quantum 
harmonic oscillator, as implied by the quasi-parabolic 
shape of the potential wells, and almost constant en- 
ergy spacing between phonon energy levels. The second, 
called the low-temperature approximation, is that only 
the lowest (zero-point) lattice vibration is excited, im- 
plying that electronic transitions do not originate from 
any of the higher phonon levels. The third, called the 
linear-coupling approximation, is that the interaction be- 
tween the defect and the lattice is the same in both the 
ground and the excited states; this implies two equally 
shaped parabolic potentials and equally spaced phonon 
energy levels in both the ground and excited states. The 
detected phonon peaks in the spectrum may be associ- 
ated with the Tfi — 1,2,... and n — 1,2,... quanta of a 
characteristic phonon mode in absorption and emission. 

The highest intensity in the phonon side band at low 
temperatures corresponds to the excitations where the 
geometry does not change, that is, the vertical absorption 




FIG. 2: (Color online) The energy (E) vs. configuration co- 
ordinate (q) diagram for the excitation process of a defect 
in the Franck-Condon approximation: E K , E c are the min- 
ima in the quasi-parabolic potential energy surfaces of the 
defect in the ground and excited states, respectively and q s , 
q e the corresponding coordinates. ZPE is the zero point en- 
ergy (indicated only for the ground state) . The energy ladders 
show the phonon energies with the phonon ground states at 
n = (ground state of the defect) and m = (excited state). 
At elevated temperatures the high-energy phonon states can 
be occupied by inducing transition A — > B (vertical absorp- 
tion, green arrow), and C — > D (vertical emission, red arrow). 
Transition A <-> C corresponds to the zero-phonon line (ZPL, 
blue double arrow) both in absorption and emission. The en- 
ergy of the Stokes shift (S) and anti-Stokes shift (AS) are also 
shown. 



and emission [l8(. This way the Stokes and anti-Stokes 
shifts are determined, which reveal the relaxation energy 
of the atoms due to the electronic excitation (see Fig. [2]). 
In the linear coupling approximation the Stokes and anti- 
Stokes shifts would have the same value but experimental 
measurements indicate a difference of about 0.05 eV (see 
Table U and related discussion) . 

Another complication in understanding the excitation 
and de-excitation processes arises from a recent measure- 
ment f]~oj ] of the highest emission intensity in the phonon 
sideband at wl.88Q eV. In this measurement, the sec- 
ond peak at 1.820 eV is cle arly visible but the third at 
1.760 eV is almost missing [10(. This measurement was 
carried out at 10 K. This is puzzling since at higher tem- 
perature (for instance, at liquid nitrogen temperature, 
77 K, where earlier experiments were conducted 0]), 
higher-energy phonon states can be occupied resulting 
in larger vertical emission energy, the opposite of what 
was found in the recent measurements [lOj. We note 
that in the latter, low-temperature (10 K) experiment, 
a laser beam was focused with a 10 cm lens, increas- 
ing the energy density by a factor of 10 4 . In fact, this 
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experiment measured the ionization of N — V~ defects 
from the 3 E excited state during the high intense exci- 
tation [lo| . Assuming that all the usual approximations 
hold for this process, an anti-Stokes shift of 0.065 eV can 
be deduced [10( which clearly contradicts the results of 
earlier experiments Q . This issue must be addressed and 
resolved in a full explanation of the process. 

State-of-the-art methods to investigate defect prop- 
erties in solids employ DFT in a supercell geometry. 
While the calculation of the ground state charge and 
spin density can be obtained accurately using the lo- 
cal density approximation (LDA) for the exchange- 
correlation functional or by other functionals that in- 
clude densi ty-g radient corrections (for example, the PBE 
functional [13}), the accurate calculation of the excita- 
tion energies presents a challenging problem due to the 
well-known self-interaction error of these methods. For 
example, the LDA value for the zero-phonon line (ZPL) 
of the N — V~ center is 1.71 eV [IH compared to the ex- 
perimental data, 1.945 eV 0- The Bethe-Salpeter equa- 
tion followed by the parameter-free GW-method for the 
quasi-particle correction of the Kohn-Sham levels [3] is 
the best tool to calculate the excitation energies, but it 
is computationally prohibitive for large supercells. Re- 
cently, it was shown that both the band gap [2(| and the 
excitation energies [2l[ are vastly improved by ap ply ing a 
screened Hartree-Fock hybrid density functional ll|, [l2j , 
referred to as the HSE06 hybrid functional. 

We employed two methods to calculate the excitation 
properties of the N — V~ center in diamond: the tra- 
ditional PBE functional and the HSE06 hybrid func- 
tional. First, the diamond primitive lattice was op- 
timized, then a simple cubic 512-atom supercell was 
constructed. Finally, we placed the negatively charged 
nitrogen-vacancy defect in the supercell, and optimized 
the structure for each given electronic configuration. We 
employed the VASP5.1 code [22j to carry out these cal- 
culations with a plane-wave basis set (using an energy 
cut-off of 420 eV) and PAW-type potentials to model the 
atomic cores [23[ [24j ■ For the optimization of the lattice 
constant we used a plane-wave cut-off energy of 840 eV 
and a 12 x 12x12 Monkhorst-Pack k-point set [25| for 
the primitive diamond lattice. For the 512-atom super- 
cell we used the T-point that provides a well-converged 
charge density. We note here that all these calculations 
use the Born-Oppenheimer approximation, so the elec- 
tron states are calculated as a function of the coordinates 
of the nuclei treated as classical particles. This approxi- 
mation is valid at low temperatures. 

The calculated lattice constant of diamond is slightly 
different in the PBE (3.567 A) and the HSE06 (3.545 A) 
approximations. The calculated band gap is very differ- 
ent in the two approaches: the PBE functional yields 
4.16 eV, while the HSE06 functional gives 5.43 eV, very 
close to the experimental value of 5.48 eV. The PBE- 
functional error is large both in absolute value (1.32 eV) 



TABLE I: The calculated vertical absorption (A — > B) and 
vertical emission energies (C — > D), and the zero-phonon line 
(ZPL) obtained with the PBE and HSE06 functionals, com- 
pared to measured values from Ref. 0]. The Stokes-shift (S) 
between the vertical absorption and the ZPL, and the anti- 
Stokes-shift (AS) between the ZPL and the vertical emission 
are also given (all values in eV). 





ZPL 


A^B 
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C -> D 


AS 


PBE 


1.706 


1.910 


0.204 


1.534 


0.172 


HSE06 


1.955 


2.213 


0.258 


1.738 


0.217 


Exp. [2] 


1.945 


2.180 


0.235 


1.760 


0.185 



and relative value (24.1%); these values are reduced by 
the HSE06 functional to 0.05 eV and 1.0%, respectively. 
We expect that the calculated excitation energies of the 
defect excited states will be also improved by using hy- 
brid functionals [H, [U HH . 

The 3 A2 ground state of the N — V~ center is ob- 
tained by spin polarized calculations both with the PBE 
and HSE06 functionals. The 3 E excited state is simply 
obtained by promoting one electron from the a\ defect 
level to the e defect level in the band gap 14 , 15| . In the 
VASP code it is possible to set the occupation numbers 
of the single particle levels, thus the 3 E excited state can 
be calculated in a self-consistent manner through such 
constrained occupation. The total energy was minimized 
for both electronic configurations as a function of the co- 
ordinates of the nuclei, which allows us to determine the 
configuration coordinates of the ground state (q g ) and 
the excited state (q e ) of the defect. The corresponding 
energy minima in the calculations are shown as E g and 
E c in Fig. O The zero-point vibration states (n = 
and m = 0) will raise these energies by a value of order 
a few tens meV, called zero-point energy (ZPE, shown 
in Fig. [2]); for example, Davies and Hamer Q deduced 
a ZPE value of «35 meV. We note that the difference 
between the ZPE values of E g and E c is expected to be 
even smaller, of order a few meV. The energy difference 
between the energy minima of E g and E e therefore gives 
a very good estimate of the ZPL (A — > C transition, see 
Fig. [5]). The transitions A — > B and C — > D are readily 
calculated by fixing the geometry at q g and q e , respec- 
tively, while varying the electronic configurations as ex- 
plained above. We note that the error associated with 
the ZPE cannot be avoided in the calculated A — > B 
and C — ► D transitions, which means that the values for 
these transitions are less accurate than for the ZPL. The 
calculated excitation energies using the PBE and HSE06 
functionals are given in Table [U 

The PBE functional gives too low value for the ZPL 
(«1.71 eV), so this gradient corrected functional does 
not improve the LDA value (1.71 eV (HI) at all. This 
is not surprising since both the LDA and the PBE func- 
tionals suffer from the self-interaction error. However, 
the HSE06 functional gives an almost perfect value, the 
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difference from experiment being smailer than 0.5%. Ap- 
parently, the HSE06 functional improves not just the 
band gap of the perfect semiconductor but the defect 
internal transition energy as well. The PBE functional 
does not improve the LDA values for the vertical absorp- 
tion energy either, which is again too low by «0.3 eV. 
The HSE06 functional yields an almost perfect value for 
the energy of the vertical absorption compared to the ex- 
perimental result (within 1.4%). The larger discrepancy 
for the A — > B transition than for the ZPL may be at- 
tributed to the intrinsic ZPE error as explained above. 
We note that the calculated Stokes shift (the relaxation 
energy), is close to the experimental result from both 
the PBE and the HSE06 functionals. The reason is that 
the self-interaction error inherent in the PBE functional 
is almost fully canceled, as the relaxation energy cor- 
responds to the energy difference between two different 
atomic configurations with the same electronic configu- 
ration. The wavefunctions and the spin density obtained 
with the HSE06 functional are somewhat more localized 
than those obtained with the PBE functional. For ex- 
ample, the integrated spin density of the 3 E state in a 
5 3 A 3 cube centered at the vacant site containing the 
three carbon atoms and the nitrogen atom is 1.61 and 
1.64 obtained in PBE and HSE06 calculations, respec- 
tively (see Fig. Q}. This is expected of the HSE06 func- 
tional which contains the Hartree-Fock exchange, giving 
more localized wave functions than the pure DFT-PBE. 

Having established the high level of accuracy of the 
HSE06 functional, we can address the issue of the verti- 
cal emission energy and the anti-Stokes shift. The value 
calculated with the PBE functional is again very low and 
not comparable to any experimental data. However, the 
HSE06 value (1.738 eV) is very close to the measurement 
of Davies and Hamer (1.760 eV) Q- Since our calcula- 
tions are valid at low temperature we assume that the low 
temperature approximation beyond the Franck-Condon 
assumption holds for this transition. We conclude that 
the anti-Stokes shift is 0.185 eV at usual experimental 
conditions at low temperatures for the N — V~ defect in 
diamond. For this defect the linear-coupling approxima- 
tion does not hold. The difference between the Stokes 
and anti-Stoke shifts is «(0.235-0.185) eV=0.050 eV in 
experiment, which compares favorably to the HSE06 
value (0.258-0.217) eV=0.041 eV. This indicates some- 
what different shape of the PES for the ground state and 
the excited state, thus different vibration modes. 

The only remaining unresolved issue is the recent ex- 
periment sug gesting a much lower anti-Stokes value of 
0.065 eV [10(. We suggest that the high energy density 
excitation in this experiment resulted in local heating 
of the sample in the area where the laser beam was fo- 
cused. The local heating of the sample will break the low 
temperature approximation and will cause a shift in the 
occupation of the phonon states from m = to m = 2. 
That may explain why the detected vertical emission en- 



ergy is larger at 10 K (at high energy density) than at 
77 K (at low energy density) excitation. Our results in- 
dicate that detailed analysis of the vibration modes of 
the 3 E excited state is important for a complete under- 
standing of the radiative emission of the N — V~ center; 
further experimental and theoretical efforts are needed in 
this direction. 
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